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Введение. Изучение магнитного состояния аустенита сталей позволило выявить механизм воздействия внешнего 
магнитного поля на сталь в процессе закалки. Предыдущие исследования установили положительное практическое 
влияние термической обработки в магнитном поле. 

Цели работы: создать компьютерную модель магнитного состояния аустенита углеродистой стали; провести вычис- 
лительные эксперименты с системой спинов при различных значениях температуры и внешнего магнитного поля. 
Материалы и методы. Использованы положения модели Изинга. Канонический ансамбль спинов моделировался 
методом Монте-Карло с использованием алгоритма Метрополиса. 

Результаты исследования. Алгоритм реализовывался при начальных параметрах, подобранных с учетом экспери- 
ментальных данных о магнитном состоянии аустенита. Изучались неоднородности этого состояния без воздействия 
магнитного поля. Получены данные о размерах ферромагнитных кластеров в аустените при различных температу- 
рах. Отмечено, что наличие внешнего магнитного поля противодействует температурному разупорядочиванию спи- 
нов. Получены данные об увеличении размеров ферромагнитных кластеров при росте напряженности магнитного 
поля. 

Обсуждение и заключения. Разработана двумерная компьютерная модель спинового состояния аустенита углероди- 
стой стали. Вычислительные эксперименты при различных параметрах модели показали, что выше температуры 
Кюри существует ближний порядок в расположении спинов. С увеличением температуры системы размеры упоря- 
доченных областей уменьшаются, а при наложении внешнего магнитного поля увеличиваются. 
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Введение. Изучается магнитное состояние аустенита сталей и влияние, которое может оказать на него 
внешнее постоянное магнитное поле. Предыдущие исследования показали, что термическая обработка в маг- 
нитном поле имеет положительное практическое значение [1-4], поэтому представляется актуальным изучение 
механизмов такого воздействия. 

К. Зинер [5] указывал на необходимость учитывать роль магнитных превращений при рассмотрении 
мартенситных переходов в сплавах на основе железа. В расчете термодинамических характеристик мартенсит- 
ного превращения он использовал два параметра — магнитную и немагнитную составляющие изменения сво- 
бодной энергии чистого железа. 

Ряд исследователей [6-8] обнаружили отклонение температурной зависимости обратной магнитной 
восприимчивости от линейного закона Кюри — Вейса в парамагнитной области некоторых сталей. При этом 
ход кривой удовлетворительно описывался функцией Ланжевена для суперпарамагнетиков. Такое изменение 
магнитных свойств аустенита сталей связано с самопроизвольными флуктуациями дальнего ферромагнитного 
порядка. 

Эксперименты [9, 10] доказали существование ближнего порядка спинов (т. н. «роев») выше темпера- 
туры Кюри с помощью магнитной дифракции нейтронов. 

Постоянное магнитное поле в процессе закалочного охлаждения воздействует на области с упорядо- 
ченным расположением спинов в аустените и на их взаимодействие. В результате инициируется зарождение 
мартенситных кристаллов на ферромагнитных кластерах. 

Материалы и методы. На данный момент существует только один точный и подробный метод экспери- 
ментальных исследований «роев» спинов — магнитная дифракция нейтронов. В данной работе магнитное со- 
стояние аустенита углеродистой стали изучено посредством вычислительного эксперимента. Применены поло- 
жения модели Изинга [11-13]. Методом Монте-Карло [14—16] смоделирован канонический ансамбль спинов, 
использование которого позволяет моделировать поведение системы спинов при постоянной температуре. Для 
получения произвольного неравномерного распределения вероятности задействован алгоритм Метрополи- 
са [17] — частный случай процедуры выборки по значимости, когда некоторые возможные выборки отбрасы- 
ваются. 

Модель Изинга используется для моделирования фазовых переходов в магнитных веществах или 
бинарных сплавах. Эта решеточная модель учитывает взаимодействия между ближайшими узлами. Спины 
представляются магнитными моментами атомов в узлах решетки, которые взаимодействуют друг с другом и 
внешним магнитным полем (при его наличии). 

Модель Изинга базируется на следующих упрощениях: 

— не учитывается кинетическая энергия узлов решетки; 

— при вычислении энергии взаимодействия спинов учитываются только ближайшие соседи; 

— предусмотрены только два возможных состояния спинов (положительное 1 или отрицательное | направление 
ВДОЛЬ ОСИ 2). 


' Канонический ансамбль — статистический ансамбль, отвечающий особой физической системе. Она обменивается энергией с окружаю- 
щей средой (термостатом), находясь с ней в тепловом равновесии, но не обменивается веществом, поскольку отделена от термостата 
непроницаемой для частиц перегородкой. Для сокращенного описания такой системы используют два параметра: число частиц М и средняя 
энергия Е [13]. 
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Как показано в [18], изучение классической двумерной модели Изинга позволяет выявить общие 
закономерности поведения магнитных систем вблизи температуры фазового перехода, даже несмотря на 
принятые упрощения. 

При использовании модели Изинга изучались такие макроскопические характеристики системы, как 
средняя энергия (Е), средняя намагниченность (М), удельная теплоемкость С и магнитная восприимчивость х. При 
вычислениях средних значений по всем конфигурациям системы учитывался период до достижения системой 
состояния релаксации, который исключался из подсчетов”. 


Результаты исследования. Опишем полную энергию системы двумерных спинов в модели Изинга с 

учетом возможного магнитного поля й: 
Е = -] 515; — ВУ 1 5. 

Здесь ./ — константа обменной связи, которая характеризует силу взаимодействия соседних спинов. № — число всех 
спинов; (1, ]) означает суммирование по всем ближайшим парам спинов. С каждым 1-м узлом решетки связано число 
5;. Оно характеризует направление магнитного момента и может принимать значения либо +1 (если спин ориенти- 
рован в положительном направлении оси 2), либо —1 (если спин ориентирован в отрицательном направлении оси 2). 

Если значение константы обменного взаимодействия больше нуля, то однонаправленное состояние двух со- 


седних спинов энергетически выгоднее, т. е. состояние с наименьшей полной энергией является ферромагнитным. В 
ином случае более предпочтительным будет состояние, когда соседние спины антипараллельны друг другу (анти- 
ферромагнитное состояние). 

Наложение внешнего магнитного поля вдоль оси 2 прибавляет или отнимает дополнительную внутреннюю 
энергию спинов соответственно их направлению по оси 2. 


В дальнейших рассуждениях о термодинамических характеристиках системы энергии ./ и й будут измерять- 
ся в единицах температуры. Это удобно при рассмотрении взаимодействия между спинами, т. к. при нагревании си- 
стемы связи между ними ослабевают. 

Найдем соотношение, связывающее удельную теплоемкость С и флуктуации энергии системы в канониче- 
ском ансамбле: 


—_ 1 2 2 
С= (Е - (Е?)). 
Намагниченность М системы вычислялась сложением всех значений 5; в данной конфигурации. Затем вы- 


числялась средняя намагниченность (М) по всем конфигурациям системы. 
Магнитная восприимчивость * при данной температуре: 
х=т (М?) — (МУ). 

Возможные конфигурации системы спинов будут определяться значениями 2^ спиновых чисел 5. Применяя 
алгоритм Метрополиса, можно сгенерировать спиновые конфигурации 5 с вероятностью и_5), а затем усреднить 
искомые величины по всем конфигурациям. Рассматривались конфигурации, отличающиеся между собой переворо- 
том одного спина. Принятие решения о перевороте того или иного спина (т. е. о принятии пробной конфигурации 5,) 
зависело от соотношения весовых функций: 

м(5р) _ бр Еб 
= о. =е т ‘ет, 
где Е; и Е, — энергии систем с конфигурациями спинов 5 и 5, соответственно. 

Спин 5; переворачивался и принималась новая конфигурация, если г > 1 или г < 1, но больше, чем слу- 
чайное число, сгенерированное с учетом равномерного распределения на отрезке от 0 до 1. В противном случае 
спин оставался неизменным. На одном шаге Монте-Карло предпринимается число попыток переворота, равное 
числу спинов системы М. 

Для двумерного случая выбраны тороидальные краевые условия: решетка представляется кольцом, в 
котором спины, находящиеся на правой границе исходной квадратной решетки, взаимодействуют со спинами, 
находящимися на левой границе. Аналогичное кольцевое взаимодействие предусмотрено для верхней и нижней 
границ. Это обеспечивает одинаковое число взаимодействий для всех спинов. Взаимодействие любого спина с 
соседями может быть рассмотрено как взаимодействие с одним спином, величина которого равна сумме вели- 
чин четырех соседних спинов (она может равняться 0, =2 или -4). В двумерном случае минимально возможное 
значение изменения энергии при опрокидывании рассматриваемого (центрального) спина составляет 4.Л. 

Алгоритм и начальные параметры реализации модели магнитного состояния аустенита. Струк- 
тура алгоритма и основные функции, разработанные при его реализации, выполнены в математическом пакете с 
учетом рекомендаций [19]. 

Модель была реализована путем последовательного выполнения перечисленных ниже операций. 














? Далее по тексту, если речь идет о конфигурациях системы, подразумеваются конфигурации после достижения системой состояния релак- 
сации. 
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1. Задание начальных условий: числа спинов системы М, константы обменного взаимодействия .Л, 
напряженности внешнего магнитного поля й, температуры 7, числа шагов метода Монте-Карло М, (аналог вре- 
мени). 

2. Выполнение функции, создающей начальную конфигурацию системы с учетом заданных №, 4, Й. 

3. Вычисление мгновенных конфигураций системы на каждый шаг №, с учетом М, ./, й и начального со- 
стояния системы. 

4. Визуализация: 

— мгновенных конфигураций системы в интересующие моменты, 

— зависимости мгновенной энергии системы от времени, 

— зависимости мгновенной намагниченности системы от времени. 

Наблюдая конфигурации спинов, можно оценить текущее магнитное состояние системы. Изучение ука- 
занных зависимостей от времени позволяет оценить необходимое время релаксации системы т для учета его в 
последующих вычислениях. 

5. Вычисление средней энергии системы (ЕЁ), удельной теплоемкости С, магнитной восприимчивости 7 
и средней намагниченности (М). 

6. Вычисление параметров, характеризующих размеры ферромагнитно-упорядоченных областей: 
средний размер, максимальный размер, минимальный размер, среднеквадратичное отклонение по размерам. 
(Вычисление идет по набору спиновых конфигураций.) 

Алгоритм реализован при начальных параметрах, подобранных с учетом экспериментальных данных о 
магнитном состоянии аустенита [8, 20]. Температура точки Кюри железоуглеродистого аустенита стали У8 со- 
ставляет -180 К, т. е. если бы удалось сохранить решетку аустенита до этой температуры, то ниже точки Кюри 
аустенит получил бы ферромагнитные свойства. Константа обменного взаимодействия двумерной модели под- 
биралась таким образом, чтобы система испытывала магнитный фазовый переход при данной температуре. 
Значение константы обменного взаимодействия для двумерного случая не совпадает с реальным значением 
этой константы для железоуглеродистого сплава. Однако в рамках двумерной модели подобраны следующие 
параметры: М = 625, Л= 0,78, й =0, №, = 2500. Благодаря этому при моделировании в интервале температур, 
включающих область точки Кюри аустенита, наблюдались эффекты, полностью аналогичные эффектам реаль- 
ной системы, испытывающей фазовый магнитный переход из парамагнитного состояния в ферромагнитное, а 
именно: 

— максимумы удельной теплоемкости и магнитной восприимчивости в точке Кюри, 
— резкое уменьшение средней намагниченности при фазовом переходе, 
— расхождение между кривой намагниченности и кривыми удельной теплоемкости и магнитной восприимчи- 


вости при переходе в ферромагнитное состояние (рис. 1). 
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Рис. 1. Поведение свойств стали У8 около ферромагнитной точки Кюри аустенита по итогам вычислительных экспериментов 


На рис. 1 видно, что намагниченность не падает до нуля в точке Кюри. Остатки «аномалий» выше точки 
Кюри неоднократно наблюдались экспериментально [6-—10, 19]. Л. Д. Ландау [12] и другие исследователи [6, 8, 9, 19] 
объясняли это флуктуациями порядка при Т> 9. Флуктуации самопроизвольной намагниченности должны особен- 
но резко сказываться вблизи точки Кюри, ибо в самой точке Кюри они стремятся к бесконечности. 
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На рис. 2 схематически изображена температурная зависимость параметра ближнего порядка с из [19]. 
б 


0 0 т 
Рис. 2. Температурная зависимость параметра олижнего магнитного порядка [19] 


Вблизи температуры Кюри с резко уменьшается, а при Т > © асимптотически стремится к нулю. Такая «не- 
резкость» температуры Кюри с точки зрения ближнего порядка предопределяет следующее явление: все аномалии 
ферромагнетиков имеют крутой максимум в точке Кюри, однако при Г> © не должны сразу идти в ноль, а лишь 
постепенно убывают в соответствии с постепенным разрушением порядка на близких расстояниях. 

Как показано в [19], энергия обменного взаимодействия часто обусловливает силы, заметно действующие 
именно на близких расстояниях. Поэтому термодинамически равновесное состояние ферромагнетика должно опре- 
деляться ближним порядком между спинами. Это частный случай общей задачи кооперативного поведения взаимо- 
действующих атомов (например, при фазовом переходе первого рода). 

Ферромагнетизм во многом квантовое явление [15]. Тем не менее даже в рамках классического статистиче- 
ского подхода теория с учетом порядка на близких расстояниях, безусловно, полезна в иллюстративном и каче- 
ственном отношении [15, 19]. 

Л. С. Стильбанс учел «магнитный порядок» на близких расстояниях в схеме Изинга [21] и теоретически 
объяснил «остатки» аномалии теплоемкости при температурах выше точки Кюри. Аномальная теплоемкость обу- 
словлена необходимостью разрушения порядка на ближних расстояниях, который частично сохраняется и выше 
точки Кюри. 

Этот метод развивал С. В. Вонсовский [22—24]. Он теоретически показал различия между ферромагнитной и 
парамагнитной точками Кюри. Итоги его работ доказывают, что парамагнитная восприимчивость вблизи точки Кю- 
ри имеет конечное значение, а не стремится к бесконечности, как это следует из теории без учета ближнего порядка. 

Такой большой интерес к упрощенной Изинговой модели объясняется тем, что обращение к ней позволяет 
разрешить принципиальные трудности, возникающие при теоретической трактовке фазовых переходов второго рода. 

Анализ областей с ближним ферромагнитным порядком в аустените с помощью вычислительного 
эксперимента. В данном разделе изучаются неоднородности магнитного состояния аустенита без внешнего воздей- 
ствия магнитным полем. 

Для визуализации мгновенной спиновой конфигурации системы в выбранные моменты спины представле- 
ны белыми и черными квадратами. Белые квадраты — спины, ориентированные в положительном направлении оси 2 
(например, под действием обменных сил и (или) приложенного внешнего магнитного поля). Черные квадраты — 
спины, ориентированные в отрицательном направлении оси 2 (рис. 3). 
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Рис. 3. Мгновенные конфигурации системы в выбранные моменты (Т= 100 К, А =0) 
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На рис. 3 представлены визуализации системы при температуре 100 К (т. е. ниже точки Кюри) без дей- 
ствия магнитного поля в различные моменты. Ниже приводится интерпретация этих конфигураций. 

— В момент времени 0 система находится в исходном состоянии, которое инициировано с помощью 
генератора случайных значений (полностью разупорядоченное состояние). 

— Переход к моменту 1 означает, что каждый спин системы (№ = 625) был рассмотрен для принятия 
решения о его перевороте согласно описанным выше условиям. 

— К третьему шагу каждый спин был трижды рассмотрен на предмет переворота. 

— Последующие моменты показывают, что система стремится к состоянию с ферромагнитным поряд- 
ком. 

— На 30-м шаге фиксируется полностью упорядоченное состояние (за исключением одного спина). 
Подобная картина (с небольшими флуктуациями 1-3 спинов) наблюдается во все последующие моменты. Всего 
шагов М, = 2500. Следовательно, система пришла в устойчивое состояние (т. е. стадия релаксация завершена 
ранее 30-го шага). 

Надо заметить, что описанная картина характерна для текущей генерации модели с заданными услови- 
ями. При формировании модели используется генератор случайных чисел, поэтому результаты пересчета будут 
отличаться в деталях, но качественно модель всегда ведет себя подобным образом при одинаковых начальных 
условиях. 

Рассмотрим графики зависимости мгновенных значений энергии (рис. 4 а) и магнитного момента 
(рис. 4 6) системы от времени. 
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Рис. 4. Зависимость мгновенных значений полной энергии системы (а) 
и магнитного момента (6) от времени 
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Из рис. 4 видно, что система быстро стремится к состоянию равновесия (процесс релаксации). После 
перехода системы в состояние равновесия отмечаются лишь небольшие флуктуации. При достижении равнове- 
сия для данной температуры по конфигурациям системы статистически вычислялись макропараметры (средняя 
энергия, намагниченность, магнитная восприимчивость, удельная теплоемкость). Для этого по соответствую- 
щим графикам оценивалось отношение времени релаксации статической системы к состоянию равновесия. От- 
секались начальные конфигурации системы в течение 400 шагов (16 % от всего времени вычисления), что с 
запасом перекрывало стадию релаксации системы, которая не превышала 150 шагов (рис. 4). 

При закалке стали У8 мартенситное превращение начинается при температуре -—500 К. Интересно ис- 
следовать поведение модели вблизи этой температуры, т. к. при закалке в магнитном поле именно здесь прояв- 
ляется действие поля, инициирующее фазовый переход. Естественно, при данной температуре устойчивым бу- 
дет парамагнитное состояние аустенита, что и показала модель. Однако в спиновом состоянии аустенита 
наблюдаются нанообласти с ближним порядком в расположении спинов [9], которые играют существенную 
роль при воздействии внешнего магнитного поля [8] (воздействие поля проявляется во влиянии на размеры и 
устойчивость таких областей, но они существуют и без внешнего поля). 

Для получения наиболее достоверных статистических данных о ферромагнитно-упорядоченных обла- 
стях в аустените моделировалась система с максимально возможным (по техническим ограничениям оператив- 
ной памяти и быстродействия процессора) числом (М№= 10000) спинов при температурах от 100 до 500 К. Дру- 
гие начальные условия остались прежними. 

На рис. 5 представлены результаты измерения упорядоченных областей в аустените. Размер определял- 
ся по количеству атомов, приходящихся на горизонтальные и вертикальные сечения упорядоченной области: [, 
— средний размер, Мах — максимальный, Мт — минимальный, со — среднеквадратичное отклонение. 
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Рис. 5. Зависимость от температуры размеров областей 
с ферромагнитным порядком в аустените 


Согласно полученным данным, флуктуации ближнего магнитного порядка в стали У8 сохраняются до 
температур около 500 К, когда при закалке стали начинается мартенситное превращение. При этом максималь- 
ные размеры таких областей в аустените при температуре 500 К могут достигать 20 атомов в сечении (-7,2 нм). 
Максимальный размер при температурах в области точки Кюри и ниже ограничен условиями модели и не пре- 
вышает 100 атомов (т. е. максимального поперечного размера моделируемой области). Наблюдаемые колебания 
на кривой максимального размера кластера связаны с использованием в модели генератора случайных чисел. 

В теории Вейса учитывается только порядок на дальних расстояниях [19], который определяется разно- 
стью количества спинов, противоположно направленных по всему домену, независимо от их детального взаим- 
ного расположения. При этом получается, что ближний порядок просто приравнивается к дальнему, т. е. числа 
ближайших соседей разной ориентации у любого спина предполагаются пропорциональными полным числам 
разнонаправленных спинов в домене. Таким образом, локальная намагниченность ближайших соседей у каждо- 
го спина равна намагниченности, рассчитанной по всему кристаллу. В реальном кристалле наблюдается другая 


картина: параллельные спины под действием обменных сил стремятся объединяться в отдельные «рои» [9, 191, 


Машино строение и машиноведение 


бл 
1 


//уезтК.допзва.га 


Вр 


Вестник Донского государственного технического университета. 2020. Т. 20, № 1. С. 51-60. 155М 1992-5980 е155М 1992-6006 
Иезший о} оп 5ие Тесйтса! Отлегзйу. 2020. Ио. 20, по. 1, рр. 51-60. 155М 1992-5980 е155М 1992-6006 





подобно атомам газа при его конденсации или образованию флуктуаций в неидеальном газе. Как видно из 
рис. 5, при низких температурах эти флуктуации очень велики. При высоких температурах более обоснованно 
приближение, не учитывающее ближний порядок. Однако и в этом случае имеются еще флуктуации параллель- 
ных спинов в малых объемах, вызываемые силами обмена и приводящие к отличной от нуля локальной намаг- 
ниченности в отсутствие результирующего магнитного момента во всем объеме (т. е. в отсутствие дальнего 
порядка). 

Таким образом, модель предсказывает существование областей с ближним магнитным порядком в 
аустените стали УЗ перед мартенситным превращением. В этих областях, с большой долей вероятности, и воз- 
никают зародыши мартенсита. Кластеры с ферромагнитным упорядочением могут накладываться на дефекты 
кристаллической решетки и области неустойчивости кристаллической решетки [4], возникающие перед фазо- 
вым переходом. В этом случае вероятность зарождения мартенсита особенно высока. 

Добавление внешнего магнитного поля в модель магнитного состояния аустенита. Наличие внеш- 
него магнитного поля противодействует температурному разупорядочиванию спинов. На рис. 6 приведены за- 
висимости параметров кластеров в аустените с ближним порядком (параметры те же, что и на рис. 5) от напря- 


женности внешнего магнитного поля для стали У8 при температуре 500 К. 
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Рис. 6. Зависимость размеров областей с ферромагнитным порядком в аустените от напряженности внешнего магнитного 
поля при температуре 500 К 


Как видно из рис. 6, с увеличением напряженности магнитного поля размеры кластеров растут вплоть 
до состояния полного магнитного упорядочивания при большой напряженности магнитного поля. Таким обра- 
зом, при закалке в магнитном поле расширяются области с ближним порядком. При этом вероятность зарожде- 
ния мартенсита выше, а число мест зарождения больше. 

Обсуждение и заключение. Разработана двумерная компьютерная модель спинового состояния аусте- 
нита углеродистой стали У8. Вычислительные эксперименты при различных параметрах модели показали, что 
выше температуры Кюри существует ближний порядок в расположении спинов. С увеличением температуры 
системы размеры упорядоченных областей уменьшаются, однако даже при температуре начала мартенситного 
превращения (500 К) возможны флуктуации ближнего порядка максимальным размером -7 нм. При наложении 
внешнего магнитного поля и усилении его напряженности увеличиваются размеры кластеров с ферромагнит- 
ным упорядочением. С помощью представленной модели в будущих работах предполагается более подробно 
исследовать влияние магнитного поля на характеристики роев спинов, их количество, время жизни и особенно- 
сти взаимодействия. 


Долгачев Ю. В. и др. К вопросу о моделировании процесса зарождения мартенсита на кластерах ферромагнитной природы 


До{еасйеу Уи. И., е!Ё . Оп то4ейпз йе татеп5йе писеапоп оп }егготазпейс сиз1егх 





Библиографический список 

1. Пустовой, В. Н. Технология бездеформационной закалки в магнитном поле тонкостенных деталей коль- 
цевой формы / В. Н. Пустовойт, Ю. В. Долгачев // Вестник Донского государственного технического университета. 
— 2011. — Т. 11, №7 (58). — С. 1064—1071. 

2. Сверхпластичность стали в температурном интервале Мд-Мн как стимул для «внутренней» магнитной 
правки / В. Н. Пустовойт, А. В. Бровер, М. Г. Магомедов, Ю. В. Долгачев // Известия высших учебных заведений. 
Северо-Кавказский регион. Технические науки. — 2006. — № 56. — С. 42—46. 

3. Разюуон, У. М. О15югНоп-Ягее Беа{ геайтеп от го4$ ш тазпейс Не!4 / У. М. Разюуон, У. У. Бовасвеу // 
Маепа[5 Регюппапсе апа Срагасептаноп. — 2019. — Уа1. 8 (2). — Р. 320-324. 
РОГ: В&рз://401.0г/10.1520/МРС20170159. 

4. Пустовойт, В. Н. Зарождение мартенсита в условиях сверхпластичности аустенита и воздействия внешне- 
го магнитного поля / В. Н. Пустовойт, Ю. В. Долгачев // Известия Волгоградского государственного технического 
университета. — 2016. — №2 (181). — С. 114-120. 

5. Гепег, С. Ппрасё оР Мазпейзт ироп МеаПагоу/ С. Гепег// ЮМ. — 1955. — Уо1.7 (5). — Р. 619-630. 
РОГ: В&рз://401.0г2/10.1007/ВЕ03377550. 

6. УогопсЫКЫт, Г. О. Апота1оиз зиреграгатаепейзта о# Фе у рБазе о 1е Ее-Сг-М№ аПоу / Г. О. УогопсыКЫт, 
Г.. №. Котазвеу, [. С. РаКаоу // Зоуте! РБуз1с$-ЗоНа Эви е. — 1975. — У]. 16 (9). — Р. 1708—1711. 

7. Ромашев, Л. Н. Изменение магнитных свойств стали вблизи мартенситной точки / Л. Н. Ромашев, 
Л. Д. Ворончихин, И. Г. Факидов // Физика металлов и металловедение. — 1973. — Т. 36, №2. — С. 291. 

8. Ризюуой, У. М. ЕепотаепенсаШу ог4еге сает$ ш аауепие аз фе агеаз оЁР тацепзИе Югтайоп / 
У. М. Ризюуой, У. У. Порвасвеу // Етегоше Маепа!5 Кезеагсв. — 2017. — Уо.6, по.2. — Р. 249-253. 
РОГ: рз://401.0ге/10.1680/Леппог.17.00042. 

9. Зроопег, 5. Эрш согге]аНопз ш поп / 5. Зроопег, В. Г. Ауегфасй // РВуз1са! Кеуеуи. — 1966. — Уо1. 142 (2). 
— Р. 291-299. РОГ: Ьфрз://401.огэ/10.1103/РБузКеу.142.291. 

10. Вохог®, В. М. Маепейс Ргорегйез оР Ме[5 апа АПоуз / В. М. Вотопй. — Сеуеапа : Атепсап Зослеу Юг 
Ме, 1959. — 349 р. 

11. 15ще, Е. Вейгае таг ТБеопе 4ез Еетотазпей$ти$ / Е. 15 // Хейзсвий Еаг РВузЦк. — 1925. — Уо|. 31 (1). 
— Р. 253—258. РОТ: В&рз://401.0го/10.1007/ВЕ02980577. 

12. Гапдам, Г.. О. Зай$Нса| Рвуз1с$. Уо1. 5 / Г. О. Гапдам, Е. М. Та и2. — ОхЮ1а : Езе\ег, 2013. — 544 р. 

13. НШ, Т. Г. Ап шбодисйоп 10 ЗайзИса! Трепподупаптсе / Т. Г.. НШ. — Мем Уокк : Соипег Согроганоп, 1986. 
— 508 р. 

14. Вадег, К. Моше Сайо Мево4$ ш ЗазНса] РНуз1с$ / К. Вшаег. — Вейт : Зрипеег Зс1епсе & Визшез$ Ме- 
Ча, 2012. — 416 р. 

15. Рузоп, Е. 7. РБазе Тгапу@оп$ ш Опапиит Зрш бузетз \Ий 150фор1с апа Мошзобор1с Пиегаснопз / 
Е.Л. узоп, Е.Н. Глеб, В. Зипоп// Лоигпа| ог ЗаНзИса РБуз1сз. — 1978. — У. 18 (4). — Р. 335-383. 
РОГ: В&рз://901.0го/10.1007/978-3-662—10018-9 12. 

16. Сгеш, М. Мегосапошса! са$ег Моше Сао знашаноп / М. Сгеи@ // Рцузса| Вемле\у ГеНегз. — 1992. — 
\Уо1. 69 (7). — Р. 1002—1005. РОГ: 6рз://901.ог/10.1103/РБузКеуГейн.69.1002. 

17. Васы, 1. ТВе плейорой$ а[еогийт / Т. Весы, Е. ЗаШхап // Сотрийпе ш Зс1епсе & Епошеепие. — 2000. — 
Уо1. 2 (1). —Р. 65—69. РОГ: В#рз://401.0г2/10.1109/5992.814660. 

18. Вонсовский, С. В. Магнетизм / С. В. Вонсовский. — Москва : Наука, 1971. — 305 с. 

19. Поршнев, С. В. Компьютерное моделирование физических процессов в пакете МАТГАВ/ 
С. В. Поршнев. — Санкт-Петербург : Лань, 2011. — 726 с. 

20. Ризюуой, У. М. Ееайгез оЁ тацепус тапзКЮоплайоп ш ее! Фитте диепсЫте ш а сопзап тазпейс ВНе]4 / 
У. М. Ризюуой, У. У. оеасвеу, Г.Р. АгеГеуа// Маепа|5 Эслепсе Еогат. — 2019. — Уо1. 946. — Р. 304-308. 
РОГ: В&рз://901.0г2/10.4028у\лу.зслепиНсле/М$Е.9 46.304. 

21. Стильбанс, Л. С. Ближний и дальний порядок в ферромагнитных телах / Л. С. Стильбанс // Журнал экс- 
периментальной и теоретической физики. — 1939. — Т.9. — С. 432. 

22. Вонсовский, С. В. О ферромагнитной и парамагнитной точках Кюри ферромагнетиков / 
С. В. Вонсовский // Доклады Академии наук СССР. — 1940. — Т. 27, № 5. — С. 550-553. 

23. Вонсовский, С. В. Ферромагнетизм как проблема упорядочения / С. В. Вонсовский // Известия Академии 
наук СССР. — Серия физическая. — 1947. — Т. 11, № 5. — С. 485-486. 

24. Вонсовский, С. В. Ферромагнетизм / С. В. Вонсовский, Я. С. Шур. — Москва ; Ленинград : ОГИЗ, 1948. 
— 816 с. 








Машино строение и машиноведение 


сл 
=) 


Бйр://уезииК.Аоп$а.га 


Вестник Донского государственного технического университета. 2020. Т. 20, № 1. С. 51-60. 155М 1992-5980 е155М 1992-6006 
Гезший о} оп Зе Тесйтса! Отоегзйу. 2020. Ио. 20, по. 1, рр. 51-60. 155М 1992-5980 е155М 1992-6006 





Поступила в редакцию 26.12.2019 
Запланирована в номер 19.02.2020 


Об авторах 

Долгачев Юрий Вячиславович, доцент кафедры «Физическое и прикладное материаловедение» 
ФГБОУ ВО «Донской государственный технический университет» (344000, РФ, г. Ростов-на-Дону, пл. Гагарина, 
1), кандидат технических наук, доцент, Кезеагсвег Г: В-2328-2016, ОВС: Вйрз://огс14.0оге/0000-0002-8558-1136, 


уии4о|(@та!.ги. 





Пустовойт Виктор Николаевич, профессор кафедры «Физическое и прикладное материаловедение» 
ФГБОУ ВО «Донской государственный технический университет» (344000, РФ, г. Ростов-на-Дону, пл. Гагари- 
на, 1), доктор технических наук, профессор, Кезеагсвег О: А-6050-2016, ОВС: В рз://огс14.0г2/0000-0001- 
6999-3520, Ирт-45(@тай.га. 








Филоненко Ирина Олеговна, студентка кафедры «Физическое и прикладное материаловедение» 
ФГБОУ ВО «Донской государственный технический университет» (344000, РФ, г. Ростов-на-Дону, пл. Гагари- 
на, 1), ОВС: № рз://огс14.0г2/0000-0002-0572-5030, 89896231497 @етай.га. 





Иванков Иван Вениаминович, студент кафедры «Физическое и прикладное материаловедение» 
ФГБОУ ВО «Донской государственный технический университет» (344000, РФ, г. Ростов-на-Дону, пл. Гагари- 
на, 1), ОВСТО: №&рз://огс14.0ге/0000-0001-6466-412Х, уапуа 1уапКоу.99 @та!.ги. 





Заявленный вклад соавторов 

Ю. В. Долгачев — формирование основной концепции, целей и задач исследования, проведение расче- 
тов, подготовка текста, формулирование выводов; В. Н. Пустовойт — научное руководство, анализ результатов 
исследований, доработка текста, корректировка выводов; И. О. Филоненко — реализация алгоритма модели в 
среде программирования, подготовка расчетов; И. В. Иванков — доработка и корректировка модели в среде 
программирования, участие в подготовке расчетов. 


Все авторы прочитали и одобрили окончательный вариант рукописи. 


